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To obtain the probability distribution of 2D crack patterns in mesoscopic regions of a disordered 
solid, the formalism of Paper I requires that a functional form associating the crack patterns (or 
states) to their formation energy be developed. The crack states are here defined by an order param- 
eter field representing both the presence and orientation of cracks at each site on a discrete square 
network. The associated Hamiltonian represents the total work required to lead an uncracked meso- 
volume into that state as averaged over the initial quenched disorder. The effect of cracks is to create 
mesovolumes having internal heterogeneity in their elastic moduli. To model the Hamiltonian, the 
effective elastic moduli corresponding to a given crack distribution are determined that includes 
crack-to-crack interactions. The interaction terms are entirely responsible for the localization tran- 
sition analyzed in Paper III. The crack-opening energies are related to these effective moduli via 
Griffith's criterion as established in Paper I. 
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INTRODUCTION 



In triaxial-stress experiments on rocks in the brittle 
regime, the onset of a macroscopic localization of defor- 
mation is usually observed around peak stress (c.f. Be- 
suelle jl| for a review) . Such departure from a macroscop- 
ically uniform deformation regime is intrinsically beyond 
the capacities of a mean-field theory, and so a specific 
model is developed here that takes the orientational na- 
ture of crack-to-crack interactions into account. 

This is the second paper in a series of three dedicated 
to exploring how the physical properties of disordered 
solids evolve as they are led to failure in a state of com- 
pression. The goal of this paper is to obtain a reasonable 
form for the Hamiltonian Ej (e, e m ) which is defined as 
the average work required to lead an intact region at 
zero deformation to the crack state denoted by j when 
the maximum applied strain is e m and where the final 
strain e is possibly different than e m due to a final un- 
loading. This Hamiltonian must be expressed in terms 
of the spatial distribution of the local order parameter 
that is the variable used to characterize the population 
of cracks in each mesovolume of a huge disordered-solid 
system. 

Most existing lattice models explore the dynamics of 
scalar order parameters either representing the break- 
down of elastic spring or beam networks under tensile 
stress |^|, or of fuse networks ||. The analogies between 
such scalar models and fracture of disordered media has 
been widely discussed M . One advantage of our approach 
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is the ability to explore interactions based on a fully ten- 
sorial description of the stress perturbations produced 
by each crack. Another is its ability to yield analyti- 
cal rather than only numerical results. Using Ej in the 
partition function established in Paper I, it is possible 
to explore the crack patterns that emerge in compressive 
settings for which isolated cracks appear in an intrinsi- 
cally stable manner no matter their size |5| , and for which 
macroscopic localization is a collective phenomenon due 
to the energetic organization of small cracks as opposed 
to an instability associated with the largest defects. In 
the present paper, we retain the leading-order effects of 
oriented-crack populations interacting in 2D. The over- 
riding importance of the long-range elastic interactions 
leaves hope that 3D generalizations would not yield qual- 
itatively different critical behavior. 



II. PRINCIPLES OF THE MODEL 

A. Order-parameter definition 

We now elaborate on the crack model introduced in 
Paper I. Each mesovolume of a huge rock system is dis- 
cretized into a square network of diamond-shaped cells of 
size A (grain sizes) and only a single crack is potentially 
present in each cell. A crack is located at the center of 
the cell and has a length d somewhere within the support 
[0, d m ] where d m is the maximum crack length (a fixed 
parameter of the system) required to satisfy d m < A. 
In the perturbative treatment of the crack interactions 
developed herein, e = (d m /A) D is taken to be a small 
number where D is the number of space-dimensions (in 
the present model, D = 2). The local order parameter 
ip(x) associated with each cell x is taken to have an am- 
plitude ip = \ip\ = (d/d m ) D and has a sign that indicates 
whether the crack is oriented at +45° or —45° relative to 
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FIG. 1: Part of the diamond network of cells that comprise a 
mesovolume. Each cell has the linear dimension A and is only 
allowed to contain one crack. The maximum length of any 
crack is d m and this length is assumed to be sufficiently small 
that (d m /A) D -C 1. The amplitude of the order parameter 
is by definition if) = \tp\ = (d/d m ) D where d is the length of 
the crack found in the cell, while the sign of <p indicates the 
orientation of the crack as shown. 



the principle-stress direction (the so-called "axial direc- 
tion"). The model is summarized in Fig. [j] 

The restriction that cracks are either at ±45° and have 
lengths less than the grain dimensions is of course a great 
simplification compared to what is found inside of real 
rocks. However, we only need to characterize the es- 
sential features of a crack population that contribute to 
localization phenomena and, to this extent, it appears 
overly complicated to model the amazing variety of crack 
geometries encountered in real rocks J^,0|. The localiza- 
tion transition involves spontaneous breaking of symme- 
tries both under translation and parity (inversion of the 
minor stress axis) as is seen from the structure of the 
shear bands emerging in the post-peak-stress regime [§. 
The essential feature of any proposed order parameter 
is that it must reflect and quantify the amount of local 
symmetry breaking and our simple model with cracks at 
either ±45° does just this. 

Furthermore, there is evidence both from acoustic- 
emissions monitoring and from direct observation af- 
ter unloading jHJ, that cracks developing prior to peak 
stress do not exceed an extent of a few grains diameters. 
This is principally because the grain contacts which break 
are much weaker than the grains and have a finite length 
so that cracks arriving in compression do so stably 
Crack coalescence is not explicitly allowed for. However, 
since several neighboring cells in a line may all contain 
cracks of the same sign, the long-range elastic effect of 
long (coalesced) cracks is effectively allowed for. Our pic- 
ture of the final shear bands experimentally observed in 
the post-peak-stress regime is that they were created by 
unstable sliding along a band weakened in the pre-peak- 



stress regime by a concentration of coherently-oriented 
cracks flsj. Our model allows small cracks to stably con- 
centrate en echelon along conjugate bands relative to the 
principal-stress direction; however, it does not model the 
final unstable rupture along a given band. 

B. Formation energy of a crack pattern 

It has been established in Paper I that to a reasonable 
approximation, the work required to form a crack state, 
as averaged over the initial disorder, separates into one 
part representing the work required to break the grain 
contacts, and a second part representing the elastic en- 
ergy stored in the cracked solid. This was expressed in 
Eqs. (26)-(28) of Paper I as 

Ej = | em :(Co-C i ):e m + ie:C J -:£. (1) 

The first term of Eq. (Q) is the energy spent in the ir- 
reversible formation of the crack state j averaged over 
quenched disorder and was obtained through an applica- 
tion of Griffith's principle. The parameter q derives from 
the quenched-disorder distribution and lies in the range 
[1/2,1] (see Sec. Ill B.2 of Paper I). The second term 
is the reversibly-stored elastic energy with Cj being the 
elastic-stiffness tensor of state j. 

Our principal task is therefore to model the way that 
cracks and collective crack states affect the overall elastic 
moduli of a mesovolume. This requires detailed knowl- 
edge of the stress (or strain) field throughout the meso- 
volume in the presence of arbitrary crack populations, 
and we treat this need using the following approxima- 
tions. First, since the cracks in the model are isolated 
one to each cell, their main effect regarding the far-field 
stress is to change the elastic moduli of their embed- 
ding cell. Such a change is modeled assuming the cracks 
to be penny-shaped ellipsoidal cavities. We ignore how 
such ellipses change shape when the applied stress is un- 
loaded/reloaded since linear elasticity alone captures the 
principle effect of how the rock becomes weaker due to 
strategic placement of cracks in cells. Since a crack oc- 
cupies a limited extent of a cell, the modification of the 
moduli is small compared to the moduli of the intact cell 
so that the resulting far-field stress field can be developed 
as a Born series. It is in the third term of this devel- 
opment that crack-to-crack interactions are first allowed 
for. Higher-order interactions (three cracks simultane- 
ously interacting and so on) are negligible to the extent 
that e = (d m /A) D can be considered small. 

III. ELASTIC ENERGY 

The goal of this section is to determine the elastic en- 
ergy E? 1 stored in a mesovolume occuping the region 
and containing the crack state j (which denotes the spa- 
tial distribution of ip(x) at all points x of f2) when a dis- 
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placement corresponding to a uniform strain tensor e' ^ The quantity we are specifically seeking to establish is 

is applied on the external surface dfl of the mesovolume. the elastic energy density E el = (r D J n \t(x) :e(x)d D x 

where we recall that I is the linear dimension of a meso- 
volume. The definitions of the strain e.y = ^(diUj+djUi) 
A. Elastic energy of a weakly heterogeneous solid an( } stress Tij = (c° + 5c ljk i)e k i give immediately the 



The effect of the crack field <f{x) is to perturb the stiff- 
ness tensor of each cell as C(x) = C° + 8c [(p(x)} where 
C° denotes the moduli of an uncracked cell (assumed 
uniform for all cells), and where 8c(x) is a small pertur- 
bation due to the possible presence of a crack as char- 
acterized by (p(x). It is established in Appendix [a| that 
the non-zero components of Sc are typically smaller than 
those of C° by a factor e = (d m /A) D <C 1. Our prob- 
lem is to resolve an elasticity boundary-value problem in 
a region containing a weakly heterogeneous stiffness 
tensor C(x). 

The displacement boundary conditions are given as 



\/x 6 dCl, u(x) = e 



(o). 



(2) 



where x denotes distance from the center of the mesovol- 
ume. Elastostatic equilibrium requires that 



CiikifydkUi + dj[8c i]k id k ui] = (3) 



throughout fl where summation over repeated indices is 
assumed both here and throughout. Due to the linearity 
of the problem, we use the elastostatic Green tensor G = 
GijXiXj for a uniform material which is a solution of 

C^djdkGimix^x^ + S^ix-x') = (4) 
Vx e 89,, G %3 {x lX ') = 0. (5) 

The components Gij(x,x') define the i th component of 
the displacement at x induced by a unit point force acting 
along the j-axis at x' . Here, Sij is the Kronecker symbol, 
and 8 D is the D-dimensional Dirac distribution. 

The solution for the displacements when no cracks are 
present is simply u(°)(a;) = e^-x throughout all of SI. 
Thus, it is a straightforward excercise to demonstrate 
that the total displacement u in the presence of the cracks 
satisfies the following integral equation 

Ui(x) = uf\x) + / G ij (x,x')d kl [Sc :jk i m dvu m ]{x')d D x' 
Jn 

(6) 

where denotes the partial derivative relative to the 
coordinate x\. Using e as the argument of a series ex- 
pansion, we write the displacements asu = u^ ' + + 
• -uW + 0(e(" +1 )), where each is 0(e m ). Collect- 
ing terms at each order of e in Eq. (||) gives the following 
recursion relation 

u [ n+1 \x) = [ G tj (x,x')d k >[Sc jk ^di>v^](xf)d D x'. 
Jn 

(7) 

The boundaries conditions used to define G guarantee 
that for all n > 0, the displacements u^ n ) are zero on the 
boundary dCl. 



following relations: 



ij ~ ^ijkl £ H + 0C ijkl£ k l 



(8) 
(9) 



i - r 

E = ^JdH / T {n - a ^:e {a) d D x (10) 

with the convention e^ 1 ) = 0. In the last expression, 
the fact that u( a ) = on the boundary for all a > 
guarantees that after integrating by parts 



r ^-a). £ {a) d D x = / T (n-a).y u (a) d D x 

n Jn 



c>n 



where we used the facts that the stress tensor is sym- 
metric and solenoidal. The n th term of the total elastic 
energy is then 



1 



E M = ± x (n). £ (0) 



(11) 



where the upper bar denotes a volume average over O. 

The first term of the elastic energy is independent of 
the tp field, and corresponds to the physically unimpor- 
tant amount of energy 



%(o) = I e (o). C (o). £ (o) 
2 



(12) 



stored in the intact state. 

For the higher orders n > 1 , is expressed by Eq. 
, and the same argument as above using the fact that 
uW = on dfl eliminates a term: 



1- 



= t £ (n-l). §c:£ (0)_ 

2 

The second term of the developement, 

= i £ (o) : jH :£ (o) 
2 



(13) 



(14) 



represents only a local dependance on 5c (and therefore 
on the crack-field) since it does not involve nested inte- 
grals over two different positions. It will be shown to rep- 
resent only the contribution of the average crack porosity 
to the stiffness of the rock. 

The third term of the development is where the desired 
crack-to-crack interactions arrive. Using the symmetry 
of Sc under the inversion of its two first or last indices, 
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and Eq. (0) to have an integral form of u^, Eq. 
transforms to 

2l D E^ = [d b u^(x)5c abcd (x)sfJd D x 



(x)d b G aj (x,x') 



xd kl 8c ]klm {x')e [ °l d D xd D x'. (15) 

This term accounts for the way that a crack present at x' 
energetically interacts with a different crack at x. This is 
the non-local interaction term that is ultimately respon- 
sible for the localization transition. The higher terms of 
the Born development can be neglected for our purposes. 



In Appendix [a[ we demonstrate that 
Sc(x) = e [Atp(x) + H\(p(x)\] 

A = | rji - 772 




B= ( -(2a -1)772 




-(2a -l)r7 2 



-m 
o 






-(1 - a)r] 3 



(19) 
(20) 

1(21) 



where (31,772,773) are positive constants expressed in Ap- 
pendix |A| in terms of the Lame parameters. 

Making the necessary contractions over the indices, we 
easily obtain the trivial (crack-independent) energy 
using Eqs. (|l2"l), ( |l6|) and (|l7|). For later convenience, 
this result is best written in matrix form as 



Elastic energy as explicit function of the crack 
field 



To establish the terms of the Born-approximated elas- 
tic energy as explicit functions of both the crack state 
(f and the imposed strain e^°\ a few definitions are first 
introduced. 

The principal axes of are along (e rac ji a i, e ax i a i) as 
denoted in Fig. [lj. Our square network of cells is rotated 
+45° from this orthonormal basis. We work here in the 
coordinates (ei,e 2 ) of the square network so that the 
applied strain takes the form 



M = lf A 7 
7 A 



(16) 



where A = e rad i a i + £ ax iai and 7 = -(e ra diai - ^axiai) are 
the imposed dilatation and shear strain. 

For convenience, we assume the intact material to be 
isotropic. Taking A + 2fi as the stress unit, where (A, jj) 
are the Lame parameters of the material, and using the 
usual tensor-to-matrix mapping of the indices (11) — > 1; 
(22) — > 2; (12) — > 3, the fourth-order stiffness tensor of 
the intact material takes the form 



1 2a - 1 
C° = [ 2a - 1 1 

l-a 



(17) 



where 



A 



A + 2/j 



(18) 



is a material-dependent constant in the range [0.5, 1]. 

The deviation <5c of this tensor due to the possible 
presence of a crack in a cell separates into an isotropic 
contribution independant of the crack's orientation, and 
into an anisotropic orientation-dependent contribution. 



25<°> = -(A, 7 ).M .(A,7) T 



M 



a 
1 - a 



Using the auxiliary field 

ip(x) = \<p(x)\ 



(22) 
(23) 

(24) 



denoting the amplitude of each crack, one similarly ob- 
tains [using Eqs. (0) and ©-©] 



(A, 7 )-M r (A, 7 ) T 



with 



Mi = -eifj 



K-2 

«3 



m 
2 

fi- 



K 2 
K 3 

2a -1 



-V2 



a)V3 



(25) 



(26) 

(27) 
(28) 



The term proportional to Tp has algebraically cancelled 
due to the symmetry of the problem under parity; inver- 
sion of the minor axis e ra di a iflips the orientation of cracks, 
and therefore changes the sign of Tp, while the energy re- 
mains necessarily unchanged. The surviving term is neg- 
ative and proportional to ip, and accounts for the soft- 
ening of the mesovolume due to the presence of cracks. 
This dependence on the total number of cracks is the only 
order-parameter dependent effect to first order. 

Last, the crack- interaction term of principal interest 
can be readily expressed from Eqs. ( jl5|) and ( |l9| ) as 



2t D E { - 2 



- f 2 (0) (0) . . 



cd^kl ■^■a.bcd-tt-ijkl faibj 

9,2,(0) (0) A R 
- Ze £ cd £ kl A abcd^i]kl9aibj 



e b cd b kl ^ n 



abcd^ijkl h>aibj 



(29) 
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where the fourth-order tensors f, g, h are functionals of ip 
and defined as 



faibj 




f d D x> 


G ai (x,x')db<pdj'<p 


(30) 


fjaibj 




( d D x' 


Gi a i}(x,x')db<pdjiip 


(31) 


haibj 




f d D x' 


G ai (x,x')dbipdjtip. 


(32) 



In the second term, the reciprocity of the Green function 
G a i{x,x') — G a i(x',x) is used as well as the notation 

G{ai} — [G a i + Gia)/2. 

The Green tensor needed here satisfies the Dirichlet 
conditions of Eq. (||) and can be obtained, in princi- 
ple, from the infinite-space Green tensor via the image 
method. However, this transforms into an infinite 
series (one term for each image), and makes the func- 
tional integrations of Paper III analytically hopeless. To 
remedy this problem, the Green function with a peri- 
odic instead of zero boundary condition is used as ersatz. 
Since u' 1 ) is only affected close to the boundaries by this 
replacement, this approximation will be considered valid 
for the evaluation of the volume inte gra l E^ . 

The double integrals of Eqs. (|30|)-(|32|) are most easily 
expressed using the 2D finite-Fourier transform 



(33) 
(34) 



F(k) = / d u xF{x)e 



where the sum over the wavevectors k is over {k = 
1-Krhijl&i \ Vi, rii 6 Z} with an upper cutoff given by 
maxrii > £/A that reflects the fact that the order- 
parameter cannot vary on scales smaller than cell sizes 
A. Since the Green function used is defined with periodic 
boundary conditions, it satisfies G a i(x, x') = G a i(x—x'). 
Its Fourier transform is easily established, and upon re- 
calling that (A + 2/x) is adopted as the stress unit, reads 



G(fc) 
k 



1 



(1 - a)k'- 



(I-akk) 



k 

W\ 



(35) 
(36) 



where I is the identity tensor. This is real and symmet- 
ric, as is G(x) itself. Since tp and i/j are real fields, one 
has tp{—k) = <p*{k) and^'(— k) = ip*(k). Using these re- 
lations, together with the identity J n d D xe tkx — i D 5 k , 
the integrals of Eqs. (^o|)-(|3^) become the following sums 



aibj 



-y 

fc#0 



\(p(k)\ 2 kjk b (5 ia - ak t k a ) (37) 



hibj = j B }_^$imkW(k)}k J k b (S la -ak l k a )(38) 



1 



pD £ 
fe#0 



Kibj = y ] \i>(k)\ 2 kjk b (5 ia - akjkg) (39) 



where 5ft denotes the real part of a complex quantity and 
where ki denotes the i th component of k = fe/||fe||. With 
the following definitions associated with the orientation 
of k 

e k = (^k) (40) 

u k = cos(2# fc ) = cos 2 6 k - sin 2 k = k\- k\ (41) 
v k = sin(26*fc) = 2 cosO k sin k = 2k\k 2 (42) 



the remaining contraction in Eq. (|29|) over the eight in- 
dices (abcdijkl) is performed. The calculation is a bit 
long but without surprise and finally produces 



(2) _ 



E 
M 2 



1 



= -(A, 7 )-M 2 .(A, 7 ) J 



a b 
b c 



(1 - a)i 2D 

where the components a, b, and c are defined 



(43) 
(44) 



fe^O 



b = ]T b k ; 

fc#0 



c = ]T c fc (45) 
fe/o 

a k = (1 - au 2 k )nl \p k \ 2 

+ 2(1 - a)u k K 1 K 2 K^fe^jfe 

+ {l-a)4\ij k \ 2 
bk = - au k v k Kin 3 ^tupk-ipkj 

+ (1 - a)v k K 2 K 3 IV'fcl 2 
c fc =(1 - av 2 k )nl \i) k \ 2 

8|) and K\ a new 



with K2,K3 defined in Eqs. 
material-dependent constant 



m - m 



(46) 



IV. SURFACE FORMATION ENERGY 



Next, we must account for the energy Ej that irre- 
versibly went into creating the cracks of a given crack 
state j at a maximum deformation e m . In Paper I, this 
contribution was obtained using Griffith's criterion as 



El = lt D Br 



: (Cq — Cj) :e Ti 



(47) 



where q derives from the quenched disorder and is 
bounded as 0.5 < q < 1. The derivation of this state- 
ment implicitly assumed that all cracks were the same 
length. In the present treatment, cracks are allowed to 
have any length in the range < d < d m . It is a straight- 
forward excercise to demonstrate that if the breaking en- 
ergies for each possible length d are all sampled from the 
same quenched-disorder distribution, then Eq. ( jl^ ) again 
holds. We forego such a demonstration. In the notation 
of the present paper we may thus state that 

= q(E^ - E el [p}) = -q(E { ^ + E^)[p] (48) 



E\ 
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where and E^> are the terms of the Born- 

development given by Eqs. (|25| ) and (£l3j) upon replacing 
the current strain parameters A and 7 by the maximum- 
achieved strain A m and j m . 



V. TEMPERATURE 

Although not required as part of the model Hamilto- 
nian, we now give an explicit in A m ,7 m approximate 
expression for the temperature by using Eq. (59) of Pa- 
per I. This temperature was derived in Paper I assuming 
only a single crack size. Unfortunately, the result does 
not easily generalize to multiple crack sizes and so we 
simply take d = d m to obtain the estimate 



£ D T{A m , ln 



(! - ffKn e l( A m,7m) 

'ln{[C/ei(A m , 7m )]9/(i-9) -1}' 



(49) 

where e\d m is how much the first-Born elastic energy in 
a mesovolume is reduced when a crack of length d m is 
introduced [c.f. Eqs. ( pBj ) and (p6|)]. The energy density 
ei(A m ,7 m ) is defined 



ei(A m ,7 m ) = i (K 2 A 2 m 



K 37„ 



(50) 



while C is a dimensionless "fracture toughness" parameter 
defined as 



c 



(A + 2fx)d n 



(51) 



There is a phase transition when (C/ei) 9 ^ 1 = 2 and T 
diverges so that all crack states become equally probable. 

We now consider whether such a phase transition is 
expected in laboratory experiments on rocks. The order 
of magnitude values T ~ 10 2 J/m 2 , d m ~ 10~ 5 m, K2 ~ 1, 
and A + 2[i ~ 10 10 Pa are appropriate for typical grains 
in rocks so that £ ~ 10~ 3 . When a rock fails in shear, 
the accumulated strain is on the order of a percent or 
two so that the maximum value of e\ of interest is also 
on the order of 10~ 4 . We thus find that at shear failure, 
C/ei ~ 10 and so we a priori expect the localization 
transition to occur prior to the temperature-divergence 
transition. This is more quantitatively demonstrated in 
Paper III. 



VI. SUMMARY 

Collecting together both the elastic energy and the sur- 
face formation energy, we obtain at last the Hamiltonian 
to be used in performing ensemble averages over crack 
states in the next paper. We write this Hamiltonian in 
the final form 



Ej(e,e m ) 
£ J (e m )b] 



£°(A, 7 ) +£ a 
-q{E™(A ri 



(A, 7 )M+£; i " t (A,7)M 
™)M+£ int (A m ,7™)M} 



where (A, 7) are the isotropic and shear strain compo- 
nents of the current strain tensor e, and (A m ,7,„) are 
similar quantities referring to the maximum achieved 
strain e m . The energy E is the trivial elastic energy 
of the uncracked state 

£°(A,7) = i{aA 2 + (l-ah 2 } 

where a is a dimensionless elastic constant in the range 
[0.5,1] defined by Eq. (|l|). 

The next term in the Born development is E av = E^ 
which depends only on the volume average if), which is 
the fraction of cracked cells in the crack state ip and is 
thus entirely independent of the spatial fluctuations of (p. 
Its dependence on the strain (A, 7) is 

£ av (A, 7 )M = -i[ K2 A 2 + K3 7 2 ]^- 

We defined e = (d m /A) D to be a small parameter where 
D = 2 is the number of space dimensions in the model, 
and d m , A, and I as respectively the linear sizes of 
the largest crack, a unit cell, and a mesovolume. The 
three coefficients m are positive dimensionless material- 
dependent constants defined by Eqs. (|27|), ( p8| ) and (^6|). 

The interaction energy E lnt = E^ involve a quadratic 
matrix operator Pk that, for each non-zero wave- vector 
fe, mixes together the Fourier modes of both ip and ip 



£ int (A, 7 )M = 

Rk = 

Ik = 
Pk = 



2{\-a)l 2D 



Y,{Rl- p k-Rk + il-P k -i k ) 



fc#0 

3?(^);»(^ fc ) 

3((^ fc ) ;3 (ijkj 

' L h M k ' 
M k A ,, 



where 5ft and 3 represent the real and imaginary part of 
a complex number. The components L k , Mk, N k de- 
pend both on the applied-strain parameters (maximum 
or actual ones), and the wave- vector k. In anticipation 
of Paper III, it is convenient to introduce 

^3 7 
oj = - 

Kl A 

as the shear-strain variable and to define the parameter 

2/i(A + n){\+2n) 



C—K 2 /Hl = l + 



\ 3 



where c > 1. The components of the matrix P k are then 

L fc (A,«) = A 2 n\{\- a ul) 
Mk(A,cu) = A 2 n\u k [(1 - a)c - av k uj] 
N k (A,u) = A 2 k\ [(l-a)c 2 

+2(1 - a)cv k uj + (1 - avl)iv 2 ] 



7 



where u k — cos(2#j.) and — sin(2#fc) are functions 
that characterize the orientation of k through its polar 
angle 8 k . 

Note that all terms contributing to the Hamiltonian 
have been written in a dimensionless form in which en- 
ergy density Ej, like stress, is measured in units of X+2/i. 



APPENDIX A: EFFECTIVE MODULI OF A 
CRACKED CELL 

A crack is modeled here as an elongated ellipse having 
a major axis of length d and a minor axes of length w in 
the limit that w/d <C 1 which corresponds to a so-called 
"penny-shaped" crack. Its long axis is by convention ori- 
ented along ei if locally ip > 0, and along §2 if <p < 0. 
The unit cell is a square whose sides are colinear with 
(ei;e2), and has a size A ;§> d since the crack is taken 
to be small. The interior of the crack is supposed to be 
much more compliant than the embedding matrix and 
all plastic deformation will be ignored; i.e., there is no 
residual stress or strain allowed for in the cracked system 
when it is unloaded to zero applied stress. 

Denoting as usual the volume average of a quantity 
with an overbar, we seek to determine the elastic-stiffness 
tensor C of a cell as defined through the relation 

Tij = CijkiSki- (Al) 

The region inside the crack is occupied by a uniform ma- 
terial of stiffness C 1 while the intact matrix surrounding 
the crack is occupied by a material of stiffness C 1 . Upon 
denoting v the volume fraction of the crack in the cell, 
we obtain directly 

Tij = (1 - v)C° m ^ kl + vCj^slt. (A2) 

Eshelby Jll| demonstrates that the strain e 1 inside an 
elliptic inclusion is uniform while Wu relates this 
strain to the strain at infinity by a tensor T 



-L rri OC 



(A3) 



With cracks considered as small inclusions in their em- 
bedding cell (i) C 1), the approximation e°° ~ e is valid 
to leading order in the above, so that 

e% = W ijkl e kl (A4) 
(l-v)W ijk i = {SikSji - vTijki) (A5) 



Using Eqs. (A3) and (A4) for the average deformation in 
and out of the inclusion, Eq. (A2) has the desired linear 
form of Eq. (Al) with an effective stiffness tensor given 

by 



J ijkl 



= c° - v (c° - r 1 1 

ijkl V, ijvnn ijmn ) 



- rnnkl 



c° 

^ijkl 



vC 



■ijmn ' 



(A6) 



This approximation is justified under the hypothesis that 
the material inside the inclusion (air) is far more compli- 
ant than the host material (solid silicate). These rela- 
tions are valid in any space dimension D . The two- 
dimensional case of interest to us here can be obtained 
from the three-dimensional Wu-Eshelby results by work- 
ing with a three-dimensional ellipsoidal inclusion having 
semi-axes of linear dimension a = dj1\ b = w/2; c = h/2 
embedded within a cell of dimension A x A x h in the 
limit that » A. In this limit, the three-dimensional 
problem becomes one in two dimensions. 

Wu expresses his tensor T in terms of a tensor S de- 
fined by Eshelby 



Ti 



1 



when i 7^ j 



2(1 2Sijij) 

Tim T1121 T1133 
T2211 T2222 T2233 
23311 T3322 T3333 

1 — Sim — £1122 —Sinz 

— 5*2211 1 — 52222 "52233 
"53311 —53322 1 — 53333 

The Eshelby [111 tensor components are defined 



(A7) 



5im — Qa I aa + RIa 

5n22 = Qb 2 I ab -RI a 
5i2i2 = §(a 2 ' - 2 



R. 



, ■ ■ b')I ab + -{I a +I b ) 



(A8) 



(A9) 
(A10) 

(All) 



with similar expressions for the remaining components 
obtained through the permutation of a, 6, c and 1, 2, 3. In 
the nota tion o f the present paper, the various parameters 
of Eqs. (|A9|)-(|AllD are defined 



Q = 



8tt(1 



3 1 - 2a r , 
rand R — 



8rr(l - cj p ) 



(A12) 



where a v = A/2 (A + fi) is the Poisson's ratio of the solid 
material (assumed isotropic), and 

du 



I a = 2nab 

I aa = 2nab 

2. 

lab 



(a 2 + u)D 



du 



{a 2 + u) 2 D 
du 



irab , 

3 J Q (a 2 + u)(b 2 + u)D 



(A13) 
(A14) 
(A15) 



with D — yj (a 2 + u)(b 2 + u) and I c = I ac = I bc = Ice 
Similar expressions are obtained for I b and hb by 
replacing a and b in the above. These elliptic integrals 
are evaluated to the leading order in the small aspect 
ratio S — b/a which gives 



I a =4n6, 7 6 = 4tt(1 - 5) 



(A16) 



in 



4tv 



-in 
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Defining parameters q and r by 
3 



q = AttQ = 



w^r dr = A * R =w^-y (A18) 



we obtain that to the leading order in 1/5 
Tnn — 
T2211 = 

?2222 = 



3 1 

4q"5 

1 - r 1 
3 ' 1 

r(H-f-r) 5 



1 



1 



r(l + f-r)£ 



(A19) 
(A20) 
(A21) 



All remaining components of T are either O(l) and there- 
fore negligible, or are unimportant for the components of 
C related to directions 1 and 2. 

To get finally the deviation Sc of the effective elastic 
moduli of the cracked cell through Eq. ( |A6|) , we note first 
that 

Att abc 2tt a 2 b tt d 2 „ tt , „ 
H= 3a = iAia = ^=^ (A22) 

where we recall that d = 2a is the crack's length, w = 2b 
its width, and 5 its aspect ratio. It is through this expres- 
sion that the small parameter e = (d m /A) 2 <C 1 enters 
the Born series. Note that ip = M = {d/d m ) 2 char- 
acterizes the extent of the crack. The third dimension 
of h = 2c goes to infinity in order to obtain the two- 
dimensional limit of this three-dimensional system. 

Replacing q and r by their expressions in terms of the 
Lame parameters A, fx, and using by convention A + 2/j 
as the stress unit, the crack-induced perturbations of the 
cell moduli are 



<5C2222 = 

5c n n = 

<5CH22 = 



-^>C 2222 T 2222 
-^Ci 122 T 22 ii 

-^C 2222 T 22 11 



TT A + 2/1 

-eV A23 

6 fi 



IT 

-—eip 



A 2 



6 fi(X + 2fx) 
tt A 

D fi 



IT A 
D fl 



<5c 2 211 = 

<5ci2i2 = — w(C < j ) 212 Ti 2 i 2 + C5 221 T 2 n 2 ) 

TT 1 H 
;CTp- 



6 2 A + fi 



<A24) 
(A25) 
(A26) 

(A27) 



with all other terms being zero except those obtained by 
the necessary symmetries under exchange of the two first 
or two last indexes. 



Using the dimensionless constants a defined in Eq. ( |Ig| ) 
and introducing the positive dimensionless coefficients r\i 



Vi = 



tt A 3 + n(\ + 2[if 
6 _ 



m = 77 



Am(A 
tt A + 2/i 
6 



2 M ) 



A 

7T A + 2/i 



12A + M' (A28) 
we obtain at last the deviation of the elastic moduli of a 

cell containing a crack with long axis oriented along ei 

(corresponding to a positive cp) 



The expression for both 
cracks is straightforward. 



772-277! -(2a -1)772 
Sc = I -(2a - 1)7/2 -m I eip 

-(l-a)r? 3 

' (A29) 

possible orientations of the 
Orienting the crack along e 2 
instead ei is equivalent to exchanging the 1 and 2 indices 
in the components of Sc which results in an exchange of 
the components Scim and <5c 2222 , all remaining compo- 
nents of 5c being unaffected by this change. Separating 
both expressions of Sc into symmetric and antisymmet- 
ric parts, and noting that Sc = trivially when ip = 
(no crack), we obtain the general expression used in Eqs. 

©-(ill). 
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